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SUMMARY 


A digital-computer simulation was made of the midcourse or ascent 
phase of a rendezvous between a ferry vehicle and a space station. The 
simulation involved a closed-loop guidance system in which both the 
relative position and relative velocity between ferry and station are 
measured (by simulated radar) and the relative -velocity corrections 
required to null the miss distance are computed and applied. The results 
are used to study the effectiveness of a particular set of guidance 
equations and to study the effects of errors in the launch conditions 
and errors in the navigation data. A number of trajectories were inves- 
tigated over a variety of initial conditions for cases in which the space 
station was in a circular orbit and also in an elliptic orbit. Trajec- 
tories are described in terms of a rotating coordinate system fixed in 
the station. 

As a result of this study the following conclusions are drawn. 
Successful rendezvous can be achieved even with launch conditions which 
are substantially less accurate than those obtained with present-day 
techniques. The average total-velocity correction required during the 
midcourse phase is directly proportional to the radar accuracy but the 
miss distance is not. Errors in the time of booster burnout or in the 
position of the ferry at booster burnout are less important than errors 
in the ferry velocity at booster burnout. The use of dead bands to 
account for errors in the navigational (radar) equipment appears to 
depend upon a compromise between the magnitude of the velocity correc- 
tions to be made and the allowable miss distance at the termination of 
the midcourse phase of the rendezvous. When approximate guidance equa- 
tions are used, there are limits on their accuracy which are dependent 
on the angular distance about the earth to the expected point of 
rendezvous. 


INTRODUCTION 


This paper is concerned with midcourse guidance of a vehicle launched 
from the earth in an attempt to rendezvous with a target in a near-earth 
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orbit. The midcourse phase is considered to start at the time of booster 
burnout (where the payload and target are several hundred miles apart) 
and to continue until the launched vehicle or payload is within 5 to 
10 miles of the target. During this time the payload is to be guided 
onto a collision course with the target. At the end of the midcourse 
phase the payload is injected into orbit in close proximity of the tar- 
get. However, this latter maneuver, generally referred to as the terminal 
phase, is treated elsewhere in the literature (for example, ref. l) and 
is not covered herein. 

During the midcourse or ascent phase, the guidance could be basically L 

inertial (ref. 2) or a closed-loop system in which both the relative posi- 
tion and the relative velocity of the payload with respect to the target 
are measured and the relative equations of motion are solved to compute 
the corrections required to null the miss distance. Reference 3 pro- 
poses this latter method for the terminal phase of rendezvous. Since 
equipment capable of measuring relative position and velocity will prob- 
ably be required for the terminal phase, such equipment is investigated 
herein for use during the midcourse phase. 

In the study of this paper, a guidance logic system based on rela- 
tive position and velocity data is developed. By simulating this guid- 
ance logic system and computing the trajectories of the target and pay- 
load, the effects of errors in launch conditions, errors in the radar 
system, dead bands in the guidance system, and so forth, are investigated 
to determine their effect on target miss distance and fuel for velocity 
corrections. In so doing it is considered that the guidance equations 
are simply a tool, and the results are generally applicable to a 
"typical" guidance system for the midcourse phase of rendezvous. 

Although "target" and "payload" are general terms in frequent use, 
it is felt that some ambiguity may result in subsequent sections. There- 
fore, in the remainder of this paper the target will be designated as a 
"space station" in orbit about the earth and the payload will be referred 
to as a "ferry vehicle. " 


SYMBOLS 


Any consistent set of units may be used. In this report. It is 
assumed that 


ge ~ 32 * 17 feet per second per second 
r e = 3 j 960 statute miles 
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1 statute mile = 5 >280 feet 
1 foot = 0. 3046 meter 

angular momentum per unit mass (Kepler's constant) 
earth gravitational constant 
height above surface of earth 
mass of ferry vehicle 

radial distance from space station to ferry vehicle 

radial distance from center of earth 

thrust 

time, measured from booster burnout 
velocity 

axes system with center fixed in space station; the Z-axis 
is perpendicular to plane of station's orbit and X- and 
Y-axes lie in plane of orbit with Y-axis always pointing 
away from center of earth 

orthogonal components of rectangular coordinate system 

angle measured about earth from some reference position 

(When station is in a circular orbit, 0 S = 0 when t = 0; 
when station is in an elliptic orbit, 0 S = 0 at perigee. ) 

angle formed by X-axis and projection of R on orbital plane 
of space station 

a determinant defined by equation (l8) 
root-mean-square error (standard deviation) 
elapsed time until rendezvous 
time until x = y = 0 or until x = y = z = 0 

time until z = 0 

angle formed by R and orbital plane of space station 
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03 


computed angular velocity 
in a circular orbit, a) 


of space station (When station is 
= 0 S = Constant. ) 


Subscripts: 
a apogee 

c computed value 

e earth 

L 

f ferry vehicle 

L launch 

o particular value of a variable at instant at which indicated 

computation is made and, thus, designates initial condition 
at each instant of time 

p perigee 

R,Rt,R& cos \|r orthogonal components of relative velocity measured 

by spherical coordinates centered in space station 

s space station 

v resultant velocity 

Primed symbols denote quantities which are in error due to radar 
inaccuracies. 

Dots over quantities denote derivatives with respect to time. 


COORDINATES AND EQUATIONS OF MOTION 


Assumptions 

In the equations of motion used in this study, the following assump- 
tions were made. The earth has an inverse -square gravitational field 
(spherically symmetric). All motions take place at an altitude where 
the effects of the air resistance can be completely neglected. Instan- 
taneous changes in the velocity can be made through the use of small 
pulses in thrust and the total fuel expenditure is sufficiently small 
to make the changes in vehicle weight negligible. 


a\-U h-i 
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Coordinates 

The motion of the orbiting space station is defined in terms of 
polar coordinates 0 S and r g with the origin fixed at the center of 
the earth. (See fig. 1.) Since the earth*s gravitational field is 
assumed to be spherically symmetric, the motion of the station is planar 
and is not dependent upon the orientation with respect to the earth of 
the orbital plane. 

The position of the ferry is defined with respect to a rectangular 
coordinate system fixed in the space station. This coordinate system, 
shown in figure 1, is constrained to move with the station at all times. 
The X- and Y-axes lie in the plane of the station* s orbit while the 
Z-axis is perpendicular to this plane. The X- and Y-axes rotate about 
the Z-axis so that the Y-axis always points away from the center of the 
earth. To accomplish this rotation the angular velocity of the X- and 
Y-axes about the Z-axis is always identical to the instantaneous angular 
velocity of the space station in its orbit about the earth. This angular 
velocity 0 S varies with time when the station is in an elliptic orbit 
and is constant (and denoted by the symbol cd) when the station is in a 
circular orbit. 

For reasons to be explained later the position of the ferry with 
respect to the rotating coordinate system of the station is also defined 
in terms of spherical coordinates. These coordinates are denoted by R, 
9, and t and their relationship to the x,y, z coordinate system is also 
shown in figure 1. 


Equations of Motion 

In this section, the equations used to define the exact motion of 
the station and ferry are given. The first equations give the position 
and velocity of the station with respect to the center of the earth and 
the second set of equations define the relative position and velocity 
of the ferry with respect to the instantaneous position of the station. 

Motions of station . - The Kepler ian equations of motion were used 
to compute the position and velocity of the space station. These equa- 
tions are, in their usual form. 



(i) 
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where C s is the constant angular momentum per unit mass of the station. 
It was found convenient to specify the conic section in terms of its 
characteristic constants and the initial (t = 0) position of the station. 
For an ellipse, the angle 0 S is defined to "be zero at perigee. In addi- 
tion, the apogee radius r s ^ a , the perigee radius r s p , and the angular 
position of the station at time zero 9 g ^ are defined. From these four 
conditions, the following relationships can be obtained: 

Hie semimajor axis 

A ’ !( r s,P * r «,a) 

the eccentricity 

r - r 
e _ x s,a r s,p 

2A 

and the angular momentum per unit mass 

C s - [ger e 2 A(l - e 2 )] l/2 (4) 

With equations (l) to (4) all necessary information on the station tra- 
jectory was obtained. 

Motions of ferry . - The equations used to define the relative motions 
of the ferry are derived in reference 4 for the general case of a sta- 
tion in an elliptic orbit. For the particular case of a station in a 
circular orbit, the equations are also given in reference 3. The general 
equations in terms of x,y,z coordinate system are 


( 2 ) 

(3) 



(5) 


(6) 
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The thrust term T/m does not appear in these equations since only 
impulsive thrust was assumed and its effect is inserted as an instan- 
taneous change in velocity. 

For the purpose of simulating the radar or other positional-rate 
measuring devices, it was found advantageous to transform equations ( 5 ) 
to spherical coordinates. To effect this change the following trans- 
formations were used: 


x = R cos i r cos 9 
y = -R cos \|r sin 0 J 
z = R sin i|r 

These transformations applied to equations ( 5 ) yield 


( 7 ) 


R - R^ - r(§ - & s )^cos ^\|r + 


2 ®e r e 

— I - S 


r s | 6 s - -^-1 - r 


r f' 


cos f sin 9 

2 r 


- (r s 0 s + 2r s 0s)cos \| |r cos 9 + — = 0 


r rt R 


i>3 


|r 9 + 2R(§ - 0 S ) - R0 s j cos ^ - 2R^(0 - 0 s )sin )|f 
+ (r s 0 s + 2r s 0 s ) sin 0 + 


r e 0, 


„ „ £ 
2 ®e r e 


s|-s - — 5 -/ “ r s 


> (8) 


cos 9=0 


2 

R\|r + 2R\|r + R^9 - 0 S ) sin f cos \|/ + (r' s © s + 2r s 9 s jsin cos 9 


A 2 Se r e 

r sl s " 3 / “ ^8 

r f J 


sin \| i cos 9=0 


J 


Likewise equation (6) becomes 


r f = (r 2 + r g ^ - 2Rr s cos \|r sin e) 


1/2 


( 9 ) 
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In the special case where the station is in a circular orbit, r c 
and S s = are constants and equations (8) reduce to the form: 


ge r e 2R ^ 
r cos \|r sin 0 + = 0 

r f 5 


R - Ri 2 - R(e - o)) 2 cos 2 \(r + |L 2 - 
£r9 + 2R(9 - a))"j cos t - 2R\j/(9 - cu)sin if + (uP - — -^ r s cos 9 =0^ (10) 

\ r f 5 / ! 


rr 

Ri(r + 2R\|f + R(9 - cousin cos \(r - [cd 2 - — ]r s s ^ n t sin 9=0 

r*-3 
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Approximate Equations Used for Guidance 

For the purpose of guiding the ferry to the station, it is desirable 
to calculate on board one of the vehicles the proper relative velocity 
required for ■ rendezvous, based on the measured instantaneous relative 
positions. Since the nonlinear differential equations used to define 
the exact motions are not particularly suitable for simple on-board com- 
putations, a simplified approximate form of these equations is used for 
guidance computations. The development of these approximate equations 
is given in both references 3 and 4 and hence is only briefly outlined 
herein. 

For the case of the station in a circular orbit, the equations of 
motion in terms of the rectangular coordinate system (eqs. (5)) may be 
linearized, and, therefore, solved in closed form, by approximating the 
gravity difference between the two vehicles. If the gravity potential 
term 


®e r e 


2 


r f 


3 



(y + r s) 2 



( 11 ) 


is expanded into a Taylor series about the origin of the rotating 
coordinate system and all terms higher than the first order are dropped, 
then the gravity difference between the two vehicles is a function of 
only one variable y and is given by the expression: 


Pi T 

e 

r f 3 
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The physical significance of this approximation is discussed in ref- 
erence 4 and it is sufficient to note here that the approximation repre- 
sents a rotating parallel gravitational field rather than the spherical 
gravitational field of equations (5). 


By expanding each of the gravitational difference terms in the 
equations of motion (eqs. (5)) and by making the approximation that r c 

and 0 S are constant (exact if the station is in a circular orbit), a 
set of linear approximate equations for nonthrusting flight is obtained: 


x - 2cqy = 0 

Q 

y + 2<nx - 3cb y = 0 
2 

z + a) z = 0 

These equations have the solutions: 


= 2 ( 2 IF * 3 y 0 )sin art: - 2 ^ cos mt + (6y 0 - 3 ^)a* + 2 §• + x 0 


( 13 ) 


= (2 -^2. - 3y 0 )cos cut + ^2- sin art: + 


- 


2 ^ 

03 


( 14 ) 


z = Z Q COS 03 1 + — sin OJt 


For any instantaneous relative position of the ferry x 0 ,y o ,z Q , it may 

be shown that the instantaneous relative velocity components required 
to achieve interception after a time axe given by the following 

expressions: 

*C _ x o( sin ^l) + sin • 1^(1 - COS 00T-|_)] 

~ u ^ — (15) 

y^C ^ X 0 * COS ^l) + y o ( k S ^ n " 3^1 COS 03T^ ^ 

— = — — (16) 
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tan 03 t^ 


( 17 ) 
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where 


A = 5^1 sin “ 8(1 “ cos <kti) 


(18) 


In tenns of the spherical coordinates R, 0 , and \|r, equations ( 15 ), 
(l6), and (17) become, respectively. 


= (k 2 + Kj sin 2 0 o - Kj^ sin 0 Q cos 0 o ^ccs 2 \|f Q + K-^ sin 2 \|f Q (19) 

^7 = ( K 1 " *2 " K 3 sin2e o + k 4 sin e o cos 9 o) sin +0 cos *0 ( 20 ) 


where 


~ = Kj sin 0 Q cos 0 Q + sin 0 O - K ^ 


( 21 ) 


K-^ = “COt CJDT^ 

Kg 

K* 


— sin gdt 1 
A i 


-*(■ 


sin (jut- - arr, cos cjut. 


"4 

"5 


. . ~ uu 1 i uji -i 

AV 1 1 

/(UT , — I — I < — UWU UJ I , I 

A^ 1 U 


= -(< 
A\ 


K|, = —(OUT, -2 + 2 COS CUT, 


— ( 1 - COS CUT-^) 


SIMULATION 


General 

Hie problem situation for the midcourse phase of the rendezvous is 
shown schematically in figure 2. Hie station is in a fixed orbit and 
the ferry is launched and guided onto an intersecting trajectory. Nav- 
igational data in the form of relative position and velocity between the 
two vehicles are obtained by radar or by some optical equipment. In 
order to simulate this physical situation, the equations of motion of 


0\-vl r H 
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the station (eqs. (l)) and the ferry (eqs. (8) or (lO)) were solved by 
using a high-speed digital computer. The computer was also used to 
generate random errors for the simulation of the ferry-borne radar and 
to simulate the ferry-borne guidance system. Each of the components 
and the information flow lines are shown in figure 3* 

To the correct relative position and velocity of the ferry are 
added random errors based on assumed characteristics which are discussed 
in the next section. The position and velocity data, modified by radar 
errors, are fed into the ferry guidance system where any corrections 
in velocity are computed. The guidance system also computes when these 
corrections should be made and triggers a switch at the desired time 
which initiates the indicated action. In the simulation the proper 
orientation and thrust control are assumed so that the velocity change 
that would have taken place is fed back into the equations of motion of 
the ferry vehicle. In the integration of these equations, the velocity 
change was made instantaneously, which implies an impulsive thrust. 

Although not included in the simulation, the attitude or orientation 
and the thrust-control channels are indicated in figure 3 with dashed 
lines. The orientation channel is placed inside the switch that activates 
the corrective thrust. This arrangement allows the vehicle to maintain 
the proper attitude for thrusting at all times. Filters can be used to 
allow only low-frequency motion and the vehicle will always be properly 
oriented at the time the corrective thrust is applied. 

At the cost of additional complication, the orientation and thrust - 
control dynamics could be included for some specific vehicle and rocket 
engine. However, it was felt that the inclusion of such items would 
mask the effects of the guidance system. 


Assumed Radar Characteristics 

In this study radar contact was assumed to exist between the sta- 
tion and the ferry. In order to arrive at a reasonable estimate for 
navigational errors, the following characteristics were ascribed to the 
system. A pulse Doppler radar utilizing a dish antenna was to be located 

on the ferry vehicle.^ - This radar was to transmit with a peak power 
output of 250 kilowatts. The return signal was given a power boost by 
means of a transponder located on the station. Measurements of the 


^Whether the radar and guidance -computer system is located on the 
ferry or on the station is only of academic interest here. If it is 
located on the station, a radio data link would be required to transmit 
the guidance information to the ferry. 
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line -of -sight range, range rate, directional angle, and angular rate 
were made at a rate of 200 samples per second and then averaged over 
l/2-second intervals. It was estimated that such a system would result 
in root-mean-square errors of 


a 9 = 


0 ^ = 30 ft in range 

O' = 3 ft /sec in range rate 
R 

= 5 X lCT 1 * radians in angle 
= 10 ^ radians/sec In angular rate 


For the simulation the distributions of these errors were assumed 
to be Gaussian. At the time of each correction, a Monte Carlo technique 
(such as that used in refs. 5 and 6) was used to pick errors for each 
of the position and velocity components,, based on the foregoing c values. 
Depending upon the sign, these errors were then aided to or subtracted 
from the true position and velocity components given by the equations of 
motion. These quantities, primed symbols in figure 3, are then fed into 
the ferry -borne guidance system. 


L 

1 

4 

7 
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GUIDANCE LOGIC SYSTEM 


General 

Equations (15) to (21) define, to a degree cf approximation, the 
relative velocity required to rendezvous as a function of the instan- 
taneous relative position of the ferry and a desired time of rendezvous. 
How these equations were utilized to effect the rendezvous and the deter- 
mination of a desired time to rendezvous is the subject of this section. 
To aid in its description a block diagram of the guidance logic system 
is shown in figure 4. Each of the component parts are described in 
following subsections. 


Calculation of Time to Intercept 

It is not unreasonable to assume that at the time of launch some 
nominal trajectory for the ferry will be specified so that near- 
rendezvous conditions will be achieved. A number of such trajectories 
are given in reference 4. It is therefore assumed in this study that 
at the time of booster burnout some approximate time to intercept is 
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denoted as ( t i) l * At subsequent time t the approximate time to 
intercept is then computed with the following equation: 

Tj_(t) » ( Ti ) l - t (22) 

As the ferry approaches the station, however, there seem to be two 
inherent disadvantages in making this approximate computation. One 
objection is that this calculation does not provide a sufficiently accu- 
rate value of upon which to base the velocity computations. The 

other objection is that it is possible for t 1 to become negative before 

rendezvous is achieved, which would result in a breakdown of the guidance 
equations. In this study it was, therefore, decided to compute over 

the final part of the trajectory from the instantaneous relative position 
and velocity of the ferry and station. The equation which was used for 
this computation was 


T l(t)| = (23) 

J optimum -R(t) 

The derivation of the expression for the desired time to Intercept is 
given in appendix A. However, as a result of the approximations which 
were made in deriving the equation, equation (23) is inaccurate for 
large values of 0 and cannot be applied over the entire range of 
Interest. In this study was computed by equation (22) until the 

value of became equal to 30° at which time equation (23) was used 

until rendezvous occurred. Physically the condition orr-^ = 30° means 

that the angle, measured about the center of the earth, between the 
instantaneous position of the station and the expected point of rendez- 
vous was 30°. 


Calculation of a) 


When the space station is in an elliptic orbit, the angular veloc- 
ity § s of the station is not a constant and, therefore, it is necessary 

to calculate a nonconstant value for o> to be used in the approximate 
guidance equations. The assumption was made that o> varies approxi- 
mately as a sinusoidal function of t. If, when 0 g = 0, od = (Up and, 

when 0 S = it, co = o^, then 


03 


^ + ^a 

— -* 

2 



cos 


(v 


% + “a 

+ — 

2 


•) 


(24) 
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When the station is in a circular orbit, a) reduces to the correct 
constant. 


Computation of Velocity Corrections 

Once the desired time to intercept has been specified, equations (15) 
to (IT) or (19) to (21) may be used to calculate the relative velocity 
components required to achieve the interception at that time. However, 
several factors must be considered when the equations are implemented 
into a guidance system (that is, guidance logic). One factor is the 
method of measuring position and velocity (navigation system). If the 
navigation system were perfect, a comparison of the computed and meas- 
ured velocity components would then yield the veiocity corrections. 
However, in practice, the radar used for navigation will have some 
inherent errors: some random and some due to alinement or bias errors. 

It is presumed that the root mean square of these errors will be known 
and the computed velocity corrections should take into account the fact 
that the measured positions and velocities may be in error to within 
some specified limits. 

A second factor involves the approximate equations used to compute 
the velocity corrections. These equations may be sufficiently accurate 
only after the two vehicles achieve a certain relative position. Thus, 
if, the computed velocity corrections have a sphere of validity, they 
should be used only in that sphere. 

A third factor involves the independence of the in-plane and out- 
of-plane velocity corrections. Since it may not be desirable or neces- 
sary to make both in-plane and out -of -plane corrections simultaneously, 
a provision for separate computations should be provided in the system. 

Taking into account these three factors, a guidance system for 
computing the velocity corrections was devised and this system is shown 
schematically in figure 4. As indicated in the figure, the radar- 
measured position of the ferry is used to compute the desired velocity 
of the ferry for rendezvous at time t-^. When the ferry is on a tra- 
jectory that will intersect the plane of the station's orbit, the time 
of this intersection is determined and is used to compute changes 

in the out-of -plane velocity. These computed velocities are compared 
with the radar -measured velocities and the differences are then deter- 
mined. The differences are multiplied by a reduction factor, to account 
for the a error of the radar, and these corrected values are then sent 
as command changes to the guidance system. The equations used for each 
of these steps and the times at which these velocity corrections are 
made are detailed in the following sections. 


ON-q 4=- h* tr 1 
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Coordinate conversion , - The computation of the desired velocity 
corrections in terms of the rectangular coordinates x,y, z or in terms 
of the spherical coordinates R, 9,i|r is one of choice. The desired 
velocities in terms of the former are given by equations ( 15 ), (l6), 
and (17). For this study it was considered to be more desirable to 
describe the motions and velocity changes in terms of the rectangular 
coordinates. Therefore, the radar -measured position of the ferry was 
first converted from spherical to rectangular coordinates by using the 
relations of equations ( 7 ). After computing the desired velocity com- 
ponents, these components were converted back to spherical form by using 
the inverse of equations ( 7 ). These inverse relations are 


= + y^~ + z^~ 


9 = tan“~L -^r 




= tan" 


* 


x2 + y2 

3 J 


( 25 ) 


These conversions are indicated in figure 4 although it should be 
emphasized that they could be eliminated by computing the desired veloc- 
ities directly in spherical coordinates. 


In -plane velocity components . - It Is desired to bring the two in- 
plane positional coordinates to zero simultaneously after some time 
interval t^. The velocities required to do this are computed by using 

equations (15) and (l6). These desired in-plane velocities were computed 
at a finite number of points along the trajectory. The choice of times 
or positions at which corrections should be made is somewhat controversial, 
since the total amount of velocity change required will be dependent upon 
how many velocity changes are made, upon the accuracy of the navigational 
equipment, and upon the line-of -sight range at the time at which the cor- 
rections are made. (See refs. 6, 7> and 8.) It was decided to make a 
correction as soon as radar acquisition could be obtained since the 
velocity required to correct a given miss distance increases as an inverse 
function of the range. It was also considered desirable to correct more 
frequently as the range decreased since the accuracy of the guidance 
equations and (possibly) that of the radar also increase as inverse func- 
tions of the range. Therefore, in this study, the initial in-plane 
velocity correction was made at a distance assumed to be that of radar 
acquisition and subsequent corrections were made each time the line-of - 
sight range was reduced by a factor of 2 (R q , R q j 2 , R q ^, R q /q, and 

so forth) until the range became less than 5 miles. This description 
is not meant to imply, however, that the last correction occurred at 
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5 miles exactly. It merely states trial the last correction occurred at 
a range greater than 5 miles and less than 10 miles. 

Out -of -plane velocity components . - When the ferry is not launched 
into the plane of the station's orbit, the motion of the ferry normal 
to the plane (herein referred to as out-of -plane motion) will be essen- 
tially sinusoidal. Computation of the desired out-of -plane velocity is 
based upon the premise that the out -of -plane position should be brought 
to zero/prior to or simultaneous with the in-plane position. The alter- 
native, that of injecting the ferry into a non-coplanar orbit and later 
making the orbits coplanar, appears to have two major disadvantages: 

(l) The guidance to match closely the in-plane velocity and position 
will be more difficult if the two vehicles are not in close proximity 
and (2) since the out -of -plane velocity correction will increase as the 
inertial velocity of the vehicle is increased, injection into orbit prior 
to making the orbits coplanar will usually be more expensive in fuel. 

The out -of -plane motion of the ferry is approximated by (see 
eqs. (14) ) 


z 


z Q cos out + — sin col - 


With the instantaneous relative position and velocity of the ferry and 
a value of ou specified, the time at which z will become zero T 2 

may be computed: 


t 2 


= i tan" 1 

(D 


ZqCC 



(26) 


There are then three possibilities: that t 2 will be greater than, 

equal to, or less than (the time when the ir -plane position will 

be zero). If at the initial or any subsequent correction point 
T 2 = T i> then no change in the out-of -plane velocity is made. If, how- 
ever, t 2 > r 1 , then a new out-of-plane velocity is computed so that 

x , y, and z will become zero simultaneously (that is, so that 
T 1 = t 2)‘ Thus, the guidance equations for the cut-of -plane velocity 
are as follows: 



o 


If 


keep 
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If z Q = 0 make z c = 0 

These three instructions are used throughout the entire rendezvous 
where the subscript o with the variable denotes the instantaneous value 
L of the variable at the time the correction is made. If, for instance, 

L at the time of radar acquisition T 2 < no correction is computed 

+ until the ferry passes through the plane of the station’s orbit (z = 0). 

I At the time the radar indicates that z = 0 , the guidance system attempts 

3 to bring z to zero. However, due to radar inaccuracies in position 

and velocity, some small error may be incurred. This error will show up 
the next time the in-plane velocity is corrected. If the newly computed 
time T 2 is then greater than the computed time t^, a second value 
of z c is computed from equation (17) which will make T 2 = t^. If, 
however, T 2 < t^, no z correction is computed until z again passes 

through zero, at which time the process is repeated. Thus, out-of -plane 
corrections are called for any time z = 0 and at any time the in -plane 
velocity is corrected, if > t-^. 

Reduction in computed velocity changes to allow for radar inaccuracy . 
As shown in figure 4, the desired values of in-plane velocity and out- 
of -plane velocity are converted to spherical coordinates and the dif- 
ference between the measured and computed values is determined. If the 
radar (excluding inaccuracies in the guidance equations) were perfectly 
accurate, these values would represent the correct velocity-component 
changes. In practice, however, the radar will be subject to errors 
which will not be removed by data smoothing and these errors may lead 
to excessive corrections in velocity. 

To compensate for such errors, so-called "dead bands" have been 
considered previously by others (for example, ref. 8). If the expected 
standard deviation of the navigational error in determining velocity <j 

is known and the computed error in velocity is less than cr v , then a 
lcTy dead band would imply that no correction should be made. If, how- 
ever, the measured velocity error were greater than the dead band, then 
(a) the total computed velocity correction could be made or (b) a veloc- 
ity correction derived from the difference between the magnitude of the 
computed velocity change and absolute magnitude of a Y could be made. 

The merits and disadvantages of such dead bands have not been completely 
established, but the principal is that of undercorrecting the computed 
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velocity change by an amount based on the mean uncertainty of navigational- 
measurement devices. In the remainder of this paper, dead bands are 
referred to as having option (a) or (b). 


In order to investigate such dead bands, it was assumed that the 
predominant errors in the radar data, after smoothing, will be due to 
"play, " hysteresis, or changes in bias and that these errors will be 
essentially statistical. It was further assumed that the root mean 
square (or standard deviation cr) of these errors was known. If a for 
each component of velocity were known and subtracted^ from the corre- 
sponding component of the corrective velocity, then the change in veloc- 
ity (with the allowance for radar errors) would be 


AV = 


■ '(M - °kf * (M - * (I 


AR0 cos \|r I - a. 


,2l 


1/2 


R0 cos if) 


However, for simplicity, it is assumed herein that the standard devia- 
tion of each velocity component is in direct proportion to the magnitude 
of the component change over the total velocity change AV C and that 

the constant of proportionality is some total expected deviation a v . 
Thus, 
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With this assumption each computed component of corrective velocity is 
reduced by the factor 


AV„ 


(|av c | > a v ) 


(27) 


The stipulation is that, if the computed change in velocity is less 
than or equal to the specified o, a zero correction is made. 
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and the total computed change in velocity is reduced by the same factor 


AV = 



(|AV C | > a v ) ( 28 ) 


As indicated in the diagram of figure 4, the difference between 
the total measured velocity and computed velocity is obtained, the fac- 
tor of equation ( 27 ) is determined, and each velocity component is 

multiplied by this factor (if cr v ^ |aV c | , a zero correction is specified 

In practice these corrections would then serve as inputs to the guidance 
system. In the simulation the velocity corrections are fed back to the 
equations of motion as instantaneous velocity changes whenever such 
changes are indicated. 


When | AV C | > a v , the foregoing procedure is essentially option (b). 

To simulate option (a), the factor of equation (27) was made unity so 
that the following relations were used in the guidance logic system: 

AV = AV C 
AV = 0 

Thus, whenever the measured error was greater than or equal to the 
predicted mean navigation error, the total computed velocity correction 
was made. If, however, the mean error was greater than or equal to the 
computed velocity change, no correction was made. 


^ | AV C | > cr v ^ 
(| AV c| = °v) 


RESULTS 


General 

In order to investigate the characteristics and effectiveness of 
the ferry guidance system employed in this study, rendezvous with a space 
station in a circular orbit and with a station in an elliptic orbit were 
attempted. At the start of each attempted rendezvous, the ferry was 
assumed to be in coasting flight (nonthrusting) after separation from 
the last stage of the boost vehicle. Nominal initial conditions were 
obtained from the studies of reference 4, Errors in the initial posi- 
tion or velocity or both were then added to simulate ferry launch errors. 
The ferry was then guided toward the station by the system described in 
the preceding section. The results and some modifications brought about 
by these results are given in the following sections. 
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Rendezvous With a Station in a Near-Earth Circular Orbit 

Typical trajectories . - With the space station in a circular orbit 
300 statute miles above the earth, rendezvous was attempted over a wide 
variety of nominal trajectories. Four of these trajectories are shown 
in figure 5* 

Presented in figures 5 (a) to 5 (d) are the variations of the range 
rate ft and the three components of the velocity in rectangular coor- 
dinates x, y> and z as a function of the range. Ticks have been used 
on each of these velocity time histories to emphasize the small dis- 
continuities due to the impulsive velocity corrections. The magnitude 
and direction of these corrections are given in a detailed breakdown in 
table I for each time history. 

Also given in these figures are the variations of the spherical 
coordinates 9 and ^ with R, the variation of the height for the 
ferry hf and station h s with distance around the earth 0, and the 
variation of the rectangular coordinates y and z with x. 

Figure 5 (a) shows a trajectory obtained when the ferry was launched 
at approximately the correct position and velocity. The initial values 
for x, y, x, and y were taken from reference 4 for a case in which 
the ferry trajectory was coplanar with the station orbit. In addition, 
initial values of 50 miles and -230 ft /sec were used for z and z, 
respectively. Any errors in the initial values of x, y, x, and y 
which result from extending the initial conditions for an in-plane 
trajectory *an out-of -plane trajectory are shorn in reference 4 to be 
relatively s'inall and therefore were not taken into account. The initial 
range from ferry to station was 457 miles. In this particular case, the 
rendezvous occurred 146 downrange from the position of the station at 
the time of ferry launch (booster burnout) and tie midcourse phase took 
38 minutes. 

The trajectory of figure 5 (b) was obtained by adding relatively 
large initial errors in position and velocity to the initial conditions 
used for figure 5 (a). The ferry was launched 8 miles too low and 
8 miles behind the desired trajectory. In addition, the ferry was 
launched with both the x and y velocity compsnents 200 ft/sec too 
low. This magnitude of error does not necessariLy represent the kind 
of accuracy that has actually been achieved in launches to date. It is 
intended to represent an extreme case. Even with such initial errors, 
rendezvous was still accomplished. 

Figures 5(c) and 5(d) show two other typical trajectories. In 
figure 5(c) the ferry approached the station from above and in fig- 
ure 5 (d) the approach is made from below the station. In both of these 
trajectories it was necessary to make a z velocity correction 
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immediately since grater than t-^. However, in figure 5(c), 

the ferry did not pass through the station's orbital plane until after 
the last correction was made (R = 5*^ miles) and the correction required 
to make the two vehicles coplanar does not appear in the time history. 

In figure 5(d), this planar intersection occurred when the ferry was 
10.6 miles from the station, and the appropriate correction was made at 
that point. 

The velocity corrections made by the simulator in performing the 
rendezvous trajectories of figure 5 are listed in table I. The table 
L lists the times in seconds at which each velocity correction was made 

1 after the initiation of the midcourse phase, the angle about the center 

4 of the earth through which the station has traveled at the time of the 

7 correction, and the range at which the velocity corrections were made, 

6 Each velocity correction is shown in terms of components in the station- 

centered rectangular coordinate system and in terms of the (station- 
centered) spherical coordinate system. The vector sum of these com- 
ponents, representing the total velocity change required at each 
correction point, is shown in the last column. The sum of the terms 
in this last column represents the total velocity correction required 
for the mission. It should be noted, however, that a la dead band was 
used in the guidance system for all the trajectories of figure 5* 

Accuracy of approximate guidance equations . - In the time histories 
of figure 5> it may be noted that no in-plane velocity correction was 
made until the range between the station and ferry was 75 to 150 miles. 
Although this line-of -sight range exceeds the upper limit of any cur- 
rently known airborne radar, the limitation was not imposed solely for 
this reason. The principal reason for the limitation is found in the 
approximate guidance equations. As pointed out in reference 4, linear- 
ization of the exact equations of motion has the practical effect of 
changing the representative spherical gravitational field of the earth 
to an approximate gravitational field. Thus, the motions computed by 
two integrations of these approximate equations becomes less accurate 
with increasing values of relative range and ® s (or cot). Furthermore, 
the guidance equations, which are the inverse of the motion equations, 
become less accurate with relative range and -0 g (or -kjdt, the distance 
to be traveled around the earth until rendezvous occurs). 

This inaccuracy can be seen in figure 6(a) which shows x, y, 

R9, and ft as computed from the guidance equations and the correct 
velocity components x, y, R9, and R required to rendezvous at each 
point on a typical trajectory (obtained from numerical solutions of the 
exact equations as described in ref. 4). The ferry trajectory lay 
entirely within the station's orbital plane. The figure shows that, 
even if the ferry is on a collision trajectory with the correct velocity, 
the guidance equations will compute that an error exists and this error 
increases with both R and orr. 
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In figure 6(b) ; the difference between the correct and computed 
values of x, y, R, and R0 is plotted as a f'uiction of out (or -fe s ). 

It may be seen that the errors are small when our < 90° but become quite 
large as gut increases beyond 90°. This effect was investigated with 
a number of trajectories where the relative ranges at the condition 
(jut = 90° were different. It was found that this effect was dependent 
primarily on out rather than on R but an Increase in R produced an 
increase in the magnitude of the error. 

When the approximate equations were used for ferry guidance at 
values of out > 90 °; it was found that large errors in the velocity 
correction were obtained. A typical result is shown in figure 6(c). 

Plotted in the figure are the velocity corrections made during an in- 
plane rendezvous where the ferry was initially on exactly the right 7 

trajectory and the correct time to rendezvous was initially specified. 6 

At the start of the trajectory, the range was 457 miles and the value 
of out was l46°. If no corrections had been made, the trajectory indi- 
cated with a solid line in figure 6(a) would have been followed. It may 
be seen that initially a very large velocity change was made followed at 
the second correction point (Rq ± 2) by an even Larger correction of the 
opposite sign. Subsequent corrections were smaller In magnitude. The 
actual trajectory (not shown) resulted in a miss of 550 feet with the 
last correction made 17 miles from the station. Even though the guidance 
system still produced a rendezvous course, the corrections were completely 
unnecessary. 

In view of these results, a limit was imposed on the guidance system 
which prevented a velocity correction from being made until the estimated 
value of outj was 90 ° or less. 

Error studies . - In order to evaluate the sensitivity of the guidance 
system to radar accuracy, a study was made of the effect of varying the 
magnitude of the assumed radar errors. This study was made by changing 
the set of values assumed for the standard deviations of the radar errors 
given in the section entitled "Assumed Radar Characteristics* " For each 
set of cf values, simulated rendezvous trajectories were computed from 
the same initial conditions a number of times, The process was performed 
with no radar error and at multiples of 1, 2, and 5 times the standard 
deviation. The sum total of the velocity corrections and the miss dis- 
tance for each trajectory were then plotted against multiples of the 
standard deviation of the radar error. The results are shown in fig- 
ure 7 as a rough envelope of the velocities required with a faired line 
through the average values. The envelope of points shows a continually 
greater spread of required velocity corrections for less accurate radar 
with a linear increase in the average value. However, the miss distance 
showed virtually no correlation or trend with radar error. In all cases 
corrections were called for at the same values of range, with the last 
velocity correction being made at 8.5 miles from the station. In this 
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study the two vehicles were coplanar, no dead band was used, and the 
velocity changes were required only to correct inaccuracies in the initial 
(launch) conditions and any errors arising from the use of the approximate 
guidance equations. 

A study was also made of varying the value of ^ initially 
specified for the guidance equations. This variation in ^ is 

equivalent to varying the position of the station at the time of booster 
burnout. For instance, in a period of 2 seconds a space station in a 
300-mile orbit would travel approximately 10 miles. The results of the 
study shown in figure 8 show virtually no trend at all for an error in 
the time of launch of up to 6 seconds. This result, at first, might seem 
surprising, but it must be realized that for the last part of the trajec- 
tory the time to intercept is computed by the guidance logic system. 

Along the first part of the flight path, a small error in the time to 
intercept does not produce a significant change in the velocity correc- 
tions. Along later parts of the flight path (<m < 30°) the time to 
intercept, as computed by the guidance logic system, is dependent only 
on the instantaneous conditions and only indirectly dependent on the 
initially specified time. Any differences in the total velocity required 
to intercept are apparently so small that they are masked by the radar 
errors. 

A study was made of the effect of varying the velocity of the ferry 
vehicle at booster burnout. Trajectories were simulated where the ferry 
was launched with an initial velocity error of ±10, ±20, and ±30 ft/sec. 

In each case, this initial error was added to the velocity components 
in rough proportion to the magnitude of the components and thus represents 
a general velocity error. The results shown in figure 9(&) show a definite 
dependence of the velocity change required upon the velocity of the ferry 
at booster burnout. 

It may be noted that the minimum of the curve of figure 9(a) occurs 
for a velocity slightly higher than the nominal value (obtained from 
ref. 4). Thus, it might appear that a slightly higher launch velocity 
is desirable. That this is not so, however, can be seen by inspection 
of figure 9(b). Figure 9(b) gives the relative closing velocity of the 
ferry and station at the time of rendezvous as a function of the error 
in launch velocity for the same trajectories of figure 9(a). It can be 
seen from this figure that a higher launch velocity also results in a 
higher closing velocity. Since this higher closing velocity must be 
removed during the terminal phase when the ferry is injected into the 
station’s orbit, there seems to be no particular advantage in a higher 
velocity at launch. In fact, when the launch velocity, the midcourse 
velocity corrections, and the velocity required to inject the ferry into 
the station’s orbit (the relative closing velocity) are added, it is 
seen that the minimum total velocity is obtained with the nominal launch 
conditions of reference 4. 
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Dead bands , - With a value of cr y of 3 ft/sec, a dead band with 
option (b) was employed in all the trajectories of the preceding sec- 
tions. Thus, no velocity correction was made if the computed total 
velocity change AV C was less than or equal to 3 ft/sec. Whenever 
AV C > 3 ft/sec, each component of the computed velocity was multiplied 
by the factor given by equation (27). Since the assumed sources of the 
errors in the radar (hysteresis, and so forth) were not functions of 
range or time, the value of cr v was kept constant throughout the entire 
trajectory. 

The effectiveness of a dead band with option (b) was evaluated by 
repeatedly launching the ferry with the same initial conditions but with 
different values of a v in the guidance logic circuit. Values of a v 
of 0, 3* 6 , 9, and 12 ft/sec were used. With each value of a v , 10 ren- 
dezvous trajectories were made. The initial conditions were the same 
as those used for the trajectory of figure 5 (a)* The average target miss 
distance and the average of the total velocity changes were computed for 
each set of 10 trajectories. The results of using a dead band with 
option (b) are shown in figure 10. The process was repeated using a 
dead band with option (a) and these results are also shown in figure 10 . 
When option (a) was used, the factor given by equation (27) was replaced 
by unity. 

The results show that the miss distance increases almost linearly 
with the magnitude of the dead band, whereas the total corrective veloc- 
ity requirements approach a limiting value. It should be noted, however, 
that a dead band of 3 ft/sec increased the miss distance by a factor of 
10 for option (b) and a factor of 5 (using the faired value) for option (a) 
over that obtained with a zero dead band. With a zero dead band the 
average miss distance was only 0.035 statute mile (185 ft) even though 
the last correction in every case was made when the ferry was 8.4 miles 
from the station. (See table l(a) which gives the relative position for 
each velocity correction. ) 


Rendezvous With a Station in a Near-Earth Elliptic Orbit 

Typical trajectories . - Rendezvous of the ferry with a station in 
an elliptic orbit (perigee 100 and apogee 500 statute miles above the 
earth) was investigated for a variety of nominal trajectories. Based 
on the launch conditions given in reference 4, rendezvous was attempted 
with the station at three positions in its orbit: 0 S = 90 °, l 80 ° (apogee), 

and 270°. Three typical trajectories, one for each of the three station 
positions at the time of rendezvous, are shown in figure 11. In addi- 
tion to the variables used in figure 5 to describe the relative motion 
of the vehicles when the station was in a circular orbit, one additional 
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variable is shown: 8> s - go. Since go is an approximate computation 

for the angular velocity (eq. (24)) of the station and is used in the 
guidance equations, a measure of the difference between 9 g and go is 

plotted for each trajectory of figure 11 as a function of the relative 
range during the rendezvous. Specific data on the magnitude and direc- 
tion of the velocity corrections applied during each trajectory are 
given in table II. 

Effect of computed value of a) . - Before arriving at the expression 
of equation (24) for computing go as a function of 0 S , a number of 
constants for go were investigated. Although 0 g varies during the 
rendezvous for a station in an elliptic orbit, rendezvous was attempted 
by using go equal to (a) the initial value of & s , (b) the expected 

final value of 0 S , and (c) the expected average value of 0 S during 
the rendezvous. Each of these approximations worked for some initial 
launch conditions and not for others. In cases where the ferry failed 
to rendezvous, several velocity corrections were made until the vehicles 
were some 20 to 30 miles apart. Before the range could be halved again, 
the range rate became positive and the nearest point of approach was some 
15 to 20 miles. If corrections had been made more often during this 
portion of' the trajectory, smaller miss distances would have been achieved. 
However, the results point up the fact that this type of an approximation 
for to is not very accurate for ellipses with an eccentricity of 0.047. 

For even smaller values of eccentricity and more frequent velocity cor- 
rections, these approximations might be sufficient. 

In all cases investigated with the station in a 100-to -500-mile 
elliptic orbit, equation (24) gave an approximate value for the angular 
velocity of the station, which was sufficiently accurate to guide the 
ferry to a rendezvous with the same accuracy as was obtained when the 
station was in a 300-mile circular orbit. It Is anticipated that more 
eccentric orbits would probably require better approximations but the 
limits have not been investigated. 

^^-i-^once lo within 200 feet of station . - Several trajectories were 
simulated with guidance corrections applied until the vehicle was between 
100 and 200 feet of the station. One typical trajectory is shown in 
figure 12 and the details of the corrections are given in table III. 

A IcJy dead band was used so that at several correction points the 
computed corrections were less than CTy and, therefore, no correction 
was made. During the last 5 miles, however, several corrections totaling 
99 ft/sec were made. In this trajectory, rendezvous occurred at about 
280 where the error in go and & s was relatively large. The station 

is accelerating as it approaches perigee and the ferry must constantly 
increase its velocity to achieve rendezvous. It may be noted that most 
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of the corrections are in the -x direction due to the acceleration of 
the station. Although R is less then one -tenth of a mile when the 
last two corrections are made, is close to 80^ so that the direc- 

tions of Ri f and -x are nearly parallel. Thus, some of these latter 
corrections have the effect of injecting the ferry into the same orbit 
as the station. The final miss distance was 16.9 feet and the injection 
velocity that would have been required at this position was 183 ft/sec. 


IMPLICATIONS OF RESULTS 

L 

1 

The foregoing results have illustrated the effectiveness of a k 

guidance technique in which the inherent errors converge to some small 7 

value or zero as the target is approached. The study is not complete 6 

in that all possible errors and facets of the guidance system were not 
investigated, but it does illustrate some of the more important gross 
effects. In addition to the finite results already given, there are 
some observations and intuitive results which merit attention here. 

The results are not especially favorable toward dead bands. How- 
ever, it should be pointed out that, if the velocity had been corrected 
every time the measured error exceeded the dead band, the results might 
have been more favorable. In the method studied, corrections were made 
only at finite intervals along the trajectory and, thus, errors were 
allowed to build up from one correction to the next. 

It is interesting to note that, when dead bands were used, most of 
the velocity corrections were applied in the same general direction; 
therefore, a ferry with a single rocket engine could have made most of 
these corrections without resorting to large angilar changes in orienta- 
tion. The velocity deficiency, usually incurred with dead bands, 
required the largest corrections to be made in tie -x directions except 
when major orbital plane changes were made. Furthermore, these correc- 
tions in the -x direction tend to inject the ferry into the station 1 s 
orbit and, hence, reduce the injection (or closing) velocity at the time 
of rendezvous. 

When the standard deviation of the radar errors was varied, the 
average miss distance was not noticeably changed but the average velocity 
requirement increased linearly with increased radar error. It would 
appear that fairly large radar errors could be tolerated during the mid- 
course phase provided velocity corrections are applied at frequent 
intervals. 

When the ferry is launched from out of the orbital plane of the 
station, the two planes may intersect before radar acquisition is 



obtained. In this case, stored or ground-transmitted data will have 
to be used to determine when and how much out-of -plane velocity should 
be applied. This correction should not represent any great difficulty. 
If, for instance, the ferry is launched out of plane but on to an 
initially parallel heading, intersection will occur 90 ° downrange from 
the launch position and a clock could be used to measure a quarter of 
the orbital period of the station. The results given herein are pos- 
sibly more applicable to errors arising from such methods since the 
error in the time interval between orbital intersection (z = 0) and the 
time at which the x component of velocity was corrected was intention- 
ally made relatively large. The digital computer could have been forced 
to pick small time intervals of integration whenever x approached zero, 
and a condition of z = 0, z = 0 could have been obtained. However, 
artifically forcing this condition would not have proved the validity 
of the guidance equations or the ability of the guidance system to bring 
z to zero prior to or simultaneous with the x and y components. 

Since the purpose of the simulation was to study gross effects, the 
angular motions of the ferry about its center of mass were not simulated 
and the thrust was applied impulsively. It is felt that the vehicle 
should be continuously oriented so that when the thrust is called for, 
it may be applied instantly without delays for reorienting the ferry. 

If a low-pass, filter were used between the guidance and orientation net- 
works, small amplitude, high-frequency excitation could be avoided. The 
assumption of impulsive thrust does not appear unreasonable since most 
of the in -plane velocity changes were on the order of 20 to 30 feet per 
second. Only the out-of -plane corrections were relatively large, being 
on the order of 100 ft/sec for every 17 miles of out-of -plane distance 
|plus any initial out-of -plane velocity) at the time of launch. However, 
these plane changes were not made at precisely z = 0 and, thus, the 
time required to make large velocity changes was probably adequatelv 
simulated. J 


CONCLUSIONS 


As a result of the foregoing studies of the midcourse phase of 
rendezvous, the following general conclusions may be drawn: 

1. The approximate guidance equations lead to large velocity errors 
in the in-plane velocity components when the earth angle until rendezvous 
( m i is more than 90 . The errors are small, however, when <m^ < 90° 

and approach zero as the angle cut^ approaches zero. This restriction 
does not apply to the out -of -plane velocity components. 
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2. Successful rendezvous are obtained even when the velocity errors 
existing at the time of booster burnout are substantially worse than can 
be expected from present-day launch accuracies- Although these errors 
are expensive in terms of velocity changes for guidance, the expense is 
not excessive. 

5- Errors in the time of booster burnout of up to 5 seconds do not 
have an appreciable effect upon the total velocity correction required. 

This result is attributed to the method of guidance which is primarily 
dependent upon the instantaneous conditions and only indirectly dependent 
upon the specified time to rendezvous. Also, for similar reasons, moderate 
errors in the launch position are relatively unimportant in the accuracy 
or fuel requirements for rendezvous. 

4 . The average total velocity correction for rendezvous is directly 
proportional to the accuracy of the radar used. However, no correlation 
is found between radar error and miss distance. 

5. When a dead band is used and corrections are made only at definite 
positions on the trajectory, the miss distance increases linearly with 
the magnitude of the dead band while the total amount of corrective veloc- 
ity approaches a limiting value. The use of dead bands appears to be 
based on a trade off between miss distance and total velocity correction. 

6. When the space station is in a near-earth elliptic orbit, ren- 
dezvous is always achieved provided a sufficiently accurate instantaneous 
value of the station angular velocity is used in the guidance equations. 

If velocity corrections are continued to within 200 feet of the station, 
most of the final correction has the effect of injecting the ferry into 
orbit. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va. , March 22, 1961. 
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APPENDIX 

COMPUTATION OF t 1 FROM RELATIVE POSITION AND 
VELOCITY OF FERRY AND STATION 


Over the final part of the trajectory the optimum time to intercept 
t-j^ was computed in the following manner. At any particular instant 

during the approach of the ferry to the space station, x Q , y Q , z Q , 
x Q , y Q , and z Q are specified (from radar or some other measurement). 

If, at that instant equations (15), (l6), and ( 17 ) are computed for a 
particular value of (along with the specified values of x Q , y Q , 

and z Q ), then the total change of velocity that must be made in order 
to intercept at time may be computed from the following equation: 


AV 




(Al) 


Now, if the process is repeated for a whole series of values of t^, 
then a variation of AV with may be plotted for that particular 

set of initial conditions (x Q , y Q , z Q , x Q , y Q , and z Q ). A typical 
plot is shown in sketch (a). 



Sketch (a) 
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It may be noted that the curve will have a number of maximums and mini- 
mums (one for each successive orbit) but the minimum of interest is the 
lowest minimum which will also usually be the first minimum. Thus, one 
method of computing the optimum time to intercept would be to generate 
such a process in the ferry vehicle and have the computer stop when the 
first minimum occurred. Another method would be to give the pilot an 
oscilloscope presentation of the plot and have him pick the minimum. In 
either case this particular value of or would be used to com- 

pute the velocity corrections that must be made to effect the rendezvous 

Because these are an infinite number of extremums, the condition 
for directly computing the extremums - namely, 



leads to a very lengthy, complicated transcendental equation. However, 
if the optimum time to intercept is "small enough" (to be defined later) 
so that 



then this condition limits the infinite number of solutions to a finite 
number of solutions, only one of which will have any physical significance. 

Such a calculation has been carried out and the condition for the 
extremums found to be 

o = Hq + a 1 (0)T 1 ) + a 2 ((UT 1 ) 2 + a^(ayr L )5 + . . . (A3) 


p P P 

^ = x o + y Q + z o 


where 
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o , J o o 

a i x o a) a) + z o a) 

Although the smallest positive roots of the cubic and quadratic expres- 
sions were computed, it was found that the simplest and best results were 
obtained from using the one unique solution obtained from the linear 
expression: 


0 = a Q + a^orr^ 


or 


T 1 X o 2 + y ° 2 + Z o 2 R o 

^opt -(x 0 x 0 + y 0 y 0 + z 0 z 0 ) 
where Rq is the line-of-sight range defined by 

% - (*o 2 * yo 2 - *o 2 f 2 

Equation ( A4) then is an expression for the optimum t ime to 
as determined from (x Q , y Q , z Q} x Q , y Q , and z Q ) or (r q 

In order to determine the accuracy of expression (a 4), 
equations (that is, from eq. (Al)) were solved for the case 
the ferry vehicle was ejected from a space station (in a 300 -mile cir- 
cular orbit) at relative velocities of 400, 600 , and 800 ft /sec at 
t = 0. The equations were solved in negative time so that the trajec- 
tories of the ferry vehicle prior to intercept were computed. Thus, 
for each value of cirrq the solutions gave the values of (x, y, z, 

x, y, and z) and (R and R). The variation of with R is 

shown in figure 13. It should be noted that oytq represents the angular 

travel of the station around the earth and as such is more significant 
than Just elapsed time. Also shown in figure 13 by dashed lines are the 
values of cut q, obtained by use of equation (A4), plotted as a function 
of R. It can be seen in figure 13 that equation (Ak) gave very good 
results out to a value of cot 1 of k0° and within the last 30 ° contin- 
ually improved as the time to intercept grew smaller. It can also be 
seen that the range of the ferry from the station was not important. 


(A4) 


intercept 
and Rq) . 


the exact 
in which 
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TABLE I. - DETAILED BREAKDOWN OF VELOCITY CORRECTIONS FROM FIGURE 5 - Continued 
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TABLE I. - DETAILED BREAKDOWN OF VELOCITY CORRECTIONS FROM FIGURE 5 - Concluded 
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TABLE II. - DETAILED BREAKDOWN OF VELOCITY CORRECTIONS FROM FIGURE 11 - Continued 
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TABUS III. - DETAILED BREAKDOWN OF VELOCITY CORRECTIONS FROM FIGURE 
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Figure 3. - Block diagram of simulation of rendezvous guidance. Components in dashed lines were 

not simulated but were assumed to exist. 








Figure k. - Block diagram of guidance system for the midcourse phase of 

rendezvous. 
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(a) Ferry launched with approximately the correct position and velocity. 
Interception occurs from below and in front of station. 


Figure 5* - Typical trajectories for rendezvous with a station in a 
300-mile circular orbit. Ferry launch occurs 50 miles out of 
plane. 
























(d) Ferry launched with small error in out-of -plane velocity. 

tion occurs from below and in front of station. 


Intercep 


Figure 5*- Concluded. 
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Exact 






(a) Magnitude of the velocity components. 

Figure 6. - The variation with position of the velocity components com- 
puted by the approximate guidance equation and the exact velocity 
components required to intercept. Variations are given in both rec- 
tangular and polar coordinates. Ferry and station are in coplanar 
orbits. 



cor, deg 


(b) Difference in magnitude between exact and computed velocity- 

components. 


Figure 6. - Continued. 

















Launch velocity error, ft/sec 


(a) Total of midcourse velocity changes as a function of launch velocity 

error. 

Figure 9. - Variation of the total velocity corrections and relative 
closing velocity with errors in the launch velocity. 
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Range, R, miles 

Figure 13* - Comparison of the exact angular distance orr^ to intercept 

(no corrections required) and the value computed from (cnR/-R). With- 
out terminal thrust control, collision would occur at relative veloc- 
ity AV. 
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